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The symplectic symmetry of eigenstates for the 0^ s in 16 O and the 0^ 3 and lowest 2 + and 4 + 
configurations of 12 C that are well-converged within the framework of the no-core shell model with 
the JISP16 realistic interaction is examined. These states are found to project at the 85 — 90% level 
onto very few symplectic representations including the most deformed configuration, which confirms 
the importance of a symplectic no-core shell model and reaffirms the relevance of the Elliott SU(3) 
model upon which the symplectic scheme is built. 
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I. INTRODUCTION 

With the development of high-precision nucleon- 
nucleon (NN) interactions derived from meson exchange 
theory and the latest advances in chiral effective field the- 
ory based on QCD that provide components of nuclear 
forces that can be matched to the underlying theory of 
quarks and gluons, ab-initio nuclear theoretical models 
play a crucial role towards a deeper understanding of 
fundamental aspects of nuclear physics. Ab-initio cal- 
culations target reproducing nuclear structure features 
while employing realistic NN (or many-nucleon) inter- 
actions and hence bridge from nuclear structure consid- 
erations to the nucleon constituent degrees of freedom 
and, in turn, to astrophysical phenomena, including nu- 
cleosynthesis and neutron stars, as well as towards a fur- 
ther exploration of 3iV nuclear forces and exotic physics 
of rare isotopes. 

The ab-initio No-Core Shell Model (NCSM) with 
modern realistic interactions yields a good description of 
the low-lying states in fcw-nuclcon systems 0] as well 
as in more complex nuclei like 12 C In addition to 
advancing our understanding of the propagation of the 
nucleon-nucleon force in nuclear matter and clustering 
phenomena [1, H|, modeling the structure of 12 C, 16 
and similar nuclei is also important for gaining a better 
understanding of other physical processes such as parity- 
violating electron scattering from light nuclei || and re- 
sults gained through neutrino studies @ as well as for 
making better predictions for capture reaction rates that 
figure prominently, for example, in the burning of He in 
massive stars Q- 

Our investigations show that the 0+ and the lowest 
2 + and 4 + states in 12 C and the 16 ground state re- 
flect the presence of an underlying symplectic sp(3,K) 



algebraic structure [17[. This is achieved through the 
projection of realistic NCSM eigenstates onto Sp(3,K)- 
symmetric basis states of the symplectic shell model that 



are free of spurious center-of-mass motion. Typically, 
eigenstates of the NCSM are reasonably well converged 
in the N max = 6 (or basis space with an effec- 

tive interaction based on the JISPI6 realistic interaction 
[1]. In particular, calculated binding energies as well 
as other observables for 12 C such as B(E2;2 1 ~ — ► 0+ ), 
B(Ml;lj" — > 0+), ground-state proton rms radii and 
the 2^ quadrupolc moment all lie reasonably close to 
the measured values. The symplectic shell model itself 
[§, [l(| is a microscopic realization of the successful Bohr- 
Mottelson collective model. It is also a multiple oscillator 
shell generalization of Elliott's SU(3) model. Hence, this 
analysis provides, for the first time, a close examination 
of symmetries in nuclei as unveiled through ab-initio cal- 
culations of the NCSM type with realistic interactions. 

The rotational nature of the ground-state band in 12 C 
nucleus has long been recognized, and early-on group- 
theoretical methods with SU(3) serving as the underpin- 
ning (sub-)symmetry were used in its description. Sym- 
plectic algebraic approaches have achieved a very good 
reproduction of low- lying energies and B(E2) values in 
light nuclei [ll|, [l2| and specifically in 12 C using phe- 
nomenological interactions (l3j or truncated symplec- 
tic basis with simplistic (semi-) microscopic interactions 
[lil . 0j5[ . Here, we establish the dominance of the sym- 
plectic Sp(3,R) symmetry in the ab-initio NCSM wave- 
functions. This in turn opens up a new and exciting 
possibility for representing significant high-M7 collective 
modes by extending the NCSM basis space beyond its 
current limits through Sp(3, R) basis states, which yields 
a dramatically smaller basis space to achieve convergence 
of higher- lying collective modes. 



II. SYMPLECTIC SHELL MODEL 

The symplectic shell model is based on the noncom- 
pact symplectic sp(3,IR) algebra that with its subalge- 
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braic structure unveils the underlying physics of a mi- 
croscopic description of collective modes in nuclei H, [l(| • 
The latter follows from the fact that the mass quadrupole 
and monopole moment operators, the many-particle ki- 
netic energy, the angular and vibrational momenta are 
all elements of the sp(3,R) D su(3) D so(3) algebraic 
structure. Hence, collective states of a nucleus with well- 
developed quadrupole and monopole vibrations as well 
as collective rotations are described naturally in terms 
of irreducible representations (irreps) of Sp(3,R). Fur- 
thermore, the elements of the sp(3,R) algebra are con- 
structed as bilinear products in the harmonic oscillator 
(HO) raising and lowering operators that in turn are ex- 
pressed through particle coordinates and linear momenta. 
This means the basis states of a Sp(3, R) irrep can be ex- 
panded in a HO (m-schemc) basis, the same basis used in 
the NCSM, thereby facilitating symmetry identification. 

The symplectic basis states are labeled (in standard 
notation [Io|) according to the reduction chain 



Sp(3,l 



D 



U(3) D 



SO(3) 
L 



(1) 



and are constructed by acting with polynomials V in the 
symplectic raising operator, A^ 2 °\ on a set of basis states 
of the symplectic bandhead, \T a ), which is a Sp(3,R) 
lowest-weight state (i?( 02 ) |r ff ) = 0, where the symplectic 
lowering operator £?( 02 ) is the adjoint of A^ 20 ^); that is, 



J k(LS)JMj 



|r CT r nP r wK (L5)JA/ J )= [p T -{A^) x |r 

(2) 

where T a = N a (A CT \i a ) labels Sp(3,R) irreps with 
(A CT fx a ) denoting a SU(3) lowest-weight state, T n = 
n(Xn^n), and T w = N u (A w /i w ). The (A„ fi n ) set gives 
the overall SU(3) symmetry of ^ coupled raising opera- 
tors in V, [X u specifies the SU(3) symmetry of the 
symplectic state, and N u = N a + n is the total number 
of oscillator quanta related to the eigenvalue, N^hH,, of a 
HO Hamiltonian that is free of spurious modes. Conse- 
quently, the symplectic basis states bring forward impor- 
tant information about the nuclear shape deformation in 
terms of the SU(3) labels, (A w /i w ), for example, (0 0), 
(A 0) and (0 fi) describe spherical, prolate and oblate 
shapes, respectively. 

The symplectic raising operator Ay®\ which is a 
SU(3) tensor with (A fi) = (2 0) character, can be ex- 
pressed as a bilinear product of the HO raising operators, 



4f = ^E< \b\xb\ 



(20) 



V2A ^s,t 



(20) 



(3) 

where the sums are over all A particles of the system. 
The first term in J3]) describes 2MI one-particle-one-hole 
(lp-lh) excitations (one particle raised by two shells) and 
the second term eliminates the spurious center-of-mass 
excitations in the construction @. For the purpose of 
comparison to NCSM results, the basis states of the \r a ) 
bandhead in (J2J) are constructed in a m-scheme basis, 

|r o .K(£ S'o)^o-^o) = 



V 



k,(L S )JoM 



10} , (4) 



where |0) is a vacuum state, Vg >il ' ) and V 
polynomials of proton (a\) and neutron (aj,) creation op- 
erators coupled to good SU(3)xSU(2) symmetry. 
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III. RESULTS AND DISCUSSIONS 

The lowest-lying eigenstates of 12 C and 16 were calcu- 
lated using the NCSM as implemented through the Many 
Fermion Dynamics (MFD) code with an effective in- 
teraction derived from the realistic JISP16 NN potential 
Q for different Ml oscillator strengths. We are particu 
larly interested in the 16 O ground state and the J 
and the lowest J = 2+(=2 ] f ) and J = 4+(=4^) states of 
the ground-state (gs) rotational band in 12 C. 



TABLE I: Probability distribution of NCSM eigenstates for 
12 C across the dominant Op-Oh and 2ftQ 2p-2h Sp(3, R) irreps, 
Ml=\h MeV. 

OhQ, 2hQ. 4hQ 6hfl Total 

J = 



Sp(3, 



NCSM 



(0 4)S = 46.26 12.58 4.76 1.24 64.84 

(1 2)5* = 1 4.80 2.02 0.92 0.38 8.12 

(1 2)5* = 1 4.72 1.99 0.91 0.37 7.99 

2hQ 2p-2h 3.46 1.02 0.35 4.83 

Total 55.78 20.05 7.61 2.34 85.78 

56.18 22.40 12.81 7.00 98.38 



Sp(3, 



NCSM 



(0 4)S = 46.80 12.41 4.55 1.19 64.95 

(1 2)5* = 1 4.84 1.77 0.78 0.30 7.69 

(1 2)5 = 1 4.69 1.72 0.76 0.30 7.47 

2hQ 2p-2h 3.28 1.04 0.38 4.70 

Total 56.33 19.18 7.13 2.17 84.81 

56.18 21.79 12.73 7.28 98.43 



Sp(3,K) (0 4)5 = 51.45 12.11 4.18 1.04 68.78 

(1 2)5 = 1 3.04 0.95 0.40 0.15 4.54 

(1 2)5 = 1 3.01 0.94 0.39 0.15 4.49 

2hQ 2p-2h 3.23 1.16 0.39 4.78 

Total 57.50 17.23 6.13 1.73 82.59 

NCSM 57.64 20.34 12.59 7.66 98.23 

For both nuclei we constructed all of the Op-Oh and 
2M7 2p-2h (2 particles raised by one shell each) symplec- 
tic bandheads and generated their Sp(3,R) irreps up to 
N m ax = 6 (6fi.f2 model space). The typical dimension 
of a symplectic irrep basis in the N max = 6 space is on 
the order of 10 2 as compared to 10 7 for the full NCSM 
m-scheme basis space. 

Analysis of overlaps of the symplectic states with the 
NCSM eigenstates for 2hVt, AMI, and 6M7 model spaces 
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FIG. 1: Ground + state probability distribution over 0M1 (blue, lowest) to 6M1 (green, highest) subspaces for the most 
dominant Op-Oh + 2hQ, 2p-2h Sp(3, R) irrep case (left) and NCSM (right) together with the leading irrep contribution (black 
diamonds), (0 4) for 12 C (a) and (0 0) for le O (b), as a function of the hQ. oscillator strength, Nmax = 6. 



TABLE II: Probability distribution of the NCSM eigenstate 
for the J = ground state in 16 O across the Op-Oh and dom- 
inant 2HQ. 2p-2h Sp(3,R) irreps, 710=15 MeV. 







Oho, 


2MI 




6m 


Total 


Sp(3,M) 


(0 0)5 = 


50.53 


15.87 


6.32 


2.30 


75.02 




2MI 2p-2h 




5.99 


2.52 


1.32 


9.83 




Total 


50.53 


21.86 


8.84 


3.62 


84.85 


NCSM 




50.53 


22.58 


14.91 


10.81 


98.83 



(N max = 2,4,6) reveals the dominance of the Op-Oh 
Sp(3, R) irreps. For the 0+ and the lowest 2+ and 4+ 
states in 12 C there are nonnegligible overlaps for only 
3 of the 13 Op-Oh Sp(3,R) irreps, namely, the leading 
(most deformed) representation with N a = 24.5 and 
(A(j fi a ) = (0 4), and carrying spin 5 — together with 
two 24.5 (12) 5 = 1 irreps with different bandhead con- 
structions for protons and neutrons. For the ground state 
of 16 there is only one possible Op-Oh Sp(3,R) irrep, 
34.5 (0 0) 5 = 0. In addition, among the 2MI 2p-2h 
Sp(3, R) irreps only a small fraction contributes signifi- 
cantly to the overlaps and it includes the most deformed 
symplectic bandhead configurations that correspond to 
oblate shapes in 12 C and prolate ones in 16 0. 

The overlaps of the most dominant symplectic states 
with the 12 C and 16 O NCSM eigenstates under consider- 
ation in the 0, 2, 4 and 6fi.il subspaces are given in TablcU 
andHIl In order to speed up the calculations, we retained 
only the largest amplitudes of the NCSM states, those 
sufficient to account for at least 98% of the norm which 
is quoted also in the table. The results show that approx- 
imately 85% of the NCSM eigenstates for 12 C ( 16 0) fall 
within a subspacc spanned by the few most significant 
Op-Oh and 2hfl 2p-2h Sp(3,R) irreps, with the 2MI 2p- 
2h Sp(3,R) irreps accounting for 5% (10%) and with the 



leading irrep, (0 4) for 12 C and (0 0) for 16 O, carrying 
close to 70% (75%) of the NCSM wavefunction. 

In addition, the projection of the NCSM wavefunc- 
tions onto the symplectic space slightly changes as one 
varies the oscillator strength Ml (see, e.g., Fig. [1] for 
the 0+. state). The overall overlaps increase towards 
smaller Ml HO frequencies and, for example, for 0+ 
it is 90% in the N max = 6 and Ml = HMeV case. 
Clearly, the largest contribution comes from the lead- 
ing, most deformed (0 4)5 = Sp(3,R) irrep for 12 C 
and (0 0)5 = for 16 0, growing to « 90% of the total 
Sp(3, R)-symmetric part for Ml =11 MeV. As expected, 
Fig. [T] also confirms that with increasing Ml the higher 
Ml excitations contribute less while the lower 0MI con- 
figurations grow in importance. 

Furthermore, the 5 = part of all three NCSM eigen- 
states for 12 C is almost entirely projected (95%) onto 
only six 5 = symplectic irreps included in Table [H 
with as much as 90% of the spin-zero NCSM states ac- 
counted solely by the leading (0 4). The 5=1 part is 
also remarkably well described by merely two Sp(3, R) ir- 
reps. Similar results are observed for the ground state of 
16 O. The outcome reveals an important property for the 
symplectic dynamical symmetry reflected within the spin 
projections of the converged NCSM states, namely, their 
Sp(3,R) symmetry and hence the geometry of the nu- 
cleon system is independent of the Ml oscillator strength 
(Fig. [J) • The symplectic symmetry is equally present in 
the spin parts of the NCSM wavefunctions for 12 C as well 
as 16 O regardless of whether the bare or the effective in- 
teractions are used. This suggests that the Lee-Suzuki 
transformation, which effectively compensates for the fi- 
nite space truncation by renormalization of the bare in- 
teraction, does not affect the Sp(3,R) symmetry struc- 
ture of the spatial wavefunctions. Hence, the symplectic 
structure detected in the present analysis for QMl model 
space is what would emerge in NSCM evaluations with a 
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FIG. 2: Projection of the S = (blue, left) [and S = 1 (red, right)] Sp(3, R) irreps onto the corresponding significant spin 
components of the NSCM wavefunctions for (a) 0^~ 3 , (b) 2+, and (c) 4^ in 12 C and (d) 0^" s in 16 0, for effective interaction for 
different HQ. oscillator strengths and bare interaction. (We present only the most significant spin values that account for more 
than 90% of a NCSM wavefunction) . 



sufficiently large model space and bare interaction. 

As N max is increased the dimension of the J = 0, 2, 
and 4 symplectic space built on the Op-Oh Sp(3,R) ir- 
reps grows very slowly compared to the NCSM space 
dimension (Fig. [3]), which remains a small fraction of 
the NCSM basis space even when the most dominant 
2hfl 2p-2h Sp(3,R) irreps are included [it is 1.29% for 
the 2hn model space, 8.7 x 1CT 2 % (4m), 8.5 x 10~ 3 % 
(6m), 9.9 x lfr 4 % (8m), 1.3 x 10- 4 % (10MI), and 
1.6 x 10 _5 % (127iQ)]. The space reduction is even more 
dramatic in the case of 16 O for the most dominant Op-Oh 
and 2MI 2p-2h configurations [0.88% for the 2M1 model 
space, 1.2xl0~ 2 % (Am), 4.4xl0~ 4 % (6hfl), 2.6xl0~ 5 % 
(8m), 2.1 x 10~ 6 % (10m), and 2.1 x 10~ 7 % (12M2)]. 
This means that a space spanned by a set of symplectic 
basis states may be computationally manageable even 
when high-M7 configurations are included. It is impor- 
tant to note that 2hil 2p-2h (2 particles raised by one 
shell each) and higher rank np-nh excitations and al- 
lowed multiples thereof can be included by building them 
into an expanded set of lowest-weight Sp(3,R) starting 
state configurations. The same "build-up" logic, 
holds because by construction these additional starting 
state configurations are also required to be lowest-weight 
Sp(3,R) states. Note that if one were to include all pos- 
sible lowest-weight np-nh starting state configurations 
(n < N max ), and allowed multiples thereof, one would 
span the entire NCSM space. 

Examination of the role of the model space truncation 
specified by N max hU, reveals that the general features of 
all outcomes are retained as the space is expanded from 
2KI to 6m (see, e.g., Fig. Hfor 0+ in 12 C). Specifically, 




S 10 ° 2 4 6 8 10 12 



N 

max 

FIG. 3: Dimension of the NCSM (blue squares) and J = 0, 2, 
and 4 Sp(3, R) (red diamonds for the 3 most significant Op-Oh 
irrep case and green circles for when all 13 Op-Oh irreps are 
included) model spaces as a function of maximum allowed 
hfl excitations, Nmax, for 12 C. 



the same three Sp(3, K) irreps in 12 C, (0 4)5 = and the 
two (1 2)5 = 1, dominate for all N max values with the 
large overlaps of the NCSM eigenstates with the lead- 
ing symplectic irreps preserved, albeit distributed out- 
ward across higher Ml excitations as the number of active 
shells increases. In this regard, it may be interesting to 
understand the importance of the latter beyond the 6M7 
model space and their role in shaping other low-lying 
states in 12 C and 16 O such as the second + . This task, 
albeit challenging, is feasible for the no-core shell model 
with symplectic Sp(3,K) extension and will be part of a 
follow-on study. 

The 0+ and 2^ states in 12 C, constructed in terms 
of the three Sp(3, R) irreps with probability ampli- 
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FIG. 4: NCSM (blue, right) and Sp(3, R) (red, left) probabil- 
ity distribution over OhQ to N max Ml subspaces for the 0^ s 
of 12 C for different model spaces, Nmax, with Ml — 15MeV. 

tudes defined by the overlaps with the NCSM wavefunc- 
tions for N max = 6 case, were also used to determine 
B(E2 : 2f -> 0+ ) transition rates. The latter, 
slightly increasing from 101% to 107% of the correspond- 
ing NCSM numbers with increasing Ml, clearly reproduce 
the NCSM results. 

The focus here has been on demonstrating the exis- 
tence of Sp(3,R) symmetry in NCSM results for 12 C, 
and therefore a possible path forward for extending the 
NCSM to a Sp-NCSM (symplectic no-core shell model) 
scheme. This will allow one to account for even higher 
Ml configurations required to realize experimentally mea- 
sured B(E2) values without an effective charge, and espe- 
cially highly deformed spatial configurations required to 
reproduce a-cluster modes in heavier nuclei. In addition, 
the results can also be interpreted as a further strong 
confirmation of Elliott's SU(3) model since the projec- 
tion of the NCSM states onto the OMl space [Fig. [TJ blue 
(right) bars] is a projection of the NCSM results onto the 
SU(3) shell model. For example, for 12 C the OMl SU(3) 
symmetry ranges from just over 40% of the NCSM 0+ 
for Ml = 11 MeV to nearly 65% for Ml =18 McV [Fig. |TJ 
blue (left) bars] with 80%-90% of this symmetry governed 
by the leading (0 4) irrep. These numbers are consistent 
with what has been shown to be a dominance of the lead- 
ing SU(3) symmetry for SU(3)-based shell- model studies 
with realistic interactions in OMl model spaces. It seems 
the simplest of Elliott's collective states can be regarded 
as a good first-order approximation in the presence of 
realistic interactions, whether the latter is restricted to 
a OMl model space or the richer multi-M] NCSM model 
spaces. 

IV. CONCLUSIONS 

In summary, we demonstrated that ab-initio NCSM 
analysis starting with the JISP16 nuclcon-nucleon inter- 



action realize a collective nuclcon motion with a clear 
symplectic symmetry structure, which moreover remains 
unaltered whether the bare or effective interactions for 
various Ml strengths are used. Specifically, NCSM wave- 
functions for the lowest 0+,, 2^~ and 4f states in 12 C 
and the ground state in 16 project at the 85-90% level 
onto very few Op-Oh and 2hfl 2p-2h spurious center-of- 
mass free symplectic irreps. While the total dimension- 
ality of the latter is only w 10~ 3 % of the NCSM space, 
they also closely reproduce the NCSM B(E2) estimates. 
The results confirm for the first time the validity of the 
Sp(3, R) approach when realistic interactions are invoked 
and when the most deformed 2MI 2p-2h symplectic ir- 
reps, which clearly improved the overlaps, arc included. 
This demonstrates the importance of the Sp(3,R) sym- 
metry, which simply matches the nuclear geometry to the 
many-nucleon dynamics, as well as reaffirm the value of 
the simpler SU(3) model upon which it is based. 

The results further suggest that a Sp-NCSM extension 
of the NCSM may be a practical scheme for achieving 
convergence to measured B(E2) values without the need 
for introducing an effective charge and even for modeling 
cluster-like phenomena as these modes can be accommo- 
dated within the general framework of the Sp(3, R) model 
if extended to large model spaces (high N max ). In addi- 
tion, the symplectic extension of the ab-initio NCSM that 
is "structured" to take advantage of massively parallel 
computing capability holds promise to allow us to model 
heavier nuclei including neutron-deficient and N ks Z nu- 
clei along the nucleosynthesis rp-path and unstable nu- 
clei currently explored in radioactive beam experiments. 
Heavy nuclei are also feasible due to the natural extension 
of the Sp(3,R) shell model to chiral-invariant pseudo- 
spin. 

In short, the NCSM with a modern realistic NN po- 
tential supports the development of collective motion in 
nuclei as can be realized within the framework of the Sp- 
NCSM and as is apparent in its Ohfl Elliott model limit. 
The Sp-NCSM is designed to model real nuclei starting 
with realistic interactions (such as the ones based on ef- 
fective field theory) , including, especially, momentum de- 
pendent forms. 
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